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We present a numerical study of polyelectrolytes electrophoresing in free solution while squeezed 
by an axisymmetric confinement force transverse to their net displacement. Hybrid multi-particle 
collision dynamics and molecular dynamics simulations with mean-field finite Debye layers show 
that even though the polyelectrolyte chains remain “free-draining”, their electrophoretic mobility 
increases with confinement in nanoconfining potential wells. The primary mechanism leading to the 
increase in mobility above the free-solution value, despite long-range hydrodynamic screening by 
counterion layers, is the orientation of polymer segments within Debye layers. The observed length- 
dependence of the electrophoretic mobility arises due to secondary effects of counterion condensation 
related to confinement compactification. 


Understanding the electrophoresis of confined polyelec¬ 
trolytes is crucial to the advancement of many separation 
techniques [T]. Methods such as translocation through 
nanopores (mi and sieving through arrays of micro¬ 
scopic posts [M9] belong to a family of techniques that 
depend on nanoengineered geometrical constraints. In 
particular, the motion of polyelectrolytes (such as DNA) 
in narrow nanochannels raises a number of fundamental 
questions that are not yet fully understood in spite of 
practical interest [iniiii]. 

It is well-known that long, electrophoresing polyelec¬ 
trolyte chains are “free-draining”: their behavior is de¬ 
scribed by local effective properties, with long-ranged 
hydrodynamic coupling mostly screened [12]. When an 
electric field E applies a force to a chain segment, it also 
applies an equal and opposite force to the diffuse layer 
of counterions over the characteristic Debye length Ad- 
Thus, the viscous forces on the surrounding fluid are ef¬ 
fectively cancelled with only a rapidly decaying residual 
hydrodynamic field beyond Ad m- This causes the effec¬ 
tive friction coefficient ^eff lo increase linearly with degree 
of polymerization N, just like the charge Qeff; the free- 
solution electrophoretic mobility /io = Qeff/Ceff is thus 
independent of chain length m- 

One way to subvert this size independence is to ap¬ 
ply a parallel mechanical force / simultaneously with the 
electric field. If the mechanical force acts only on the 
monomers and not on the counterions then the net force 
on the fluid does not cancel, long-ranged flows are pos¬ 
sible and the chain is no longer free-draining [15]. This 
concept has been applied to many situations [12] includ¬ 
ing tethered polyelectrolytes m and end-labeled free- 
solution electrophoresis m- One oft-given example is 
the electrophoretic motion of a deformed polyelectrolyte 
through a nanofluidic channel [MEn]: It is argued that 
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FIG. 1. A negatively charged polyelectrolyte chain under¬ 
going electrophoretic motion under the action of an electric 
field E. Point-like MPCD fluid particles are assigned a charge 
specified by the Debye-Hiickel approximation and feel an elec¬ 
tric force in the opposite direction to the monomer. A cylin¬ 
drical harmonic potential radially confines the chain, creating 
an effective tube of radius i^eff = yJk^T/k. 


confinement increases the frequency of collisions between 
the monomers and the channel walls, thereby inducing 
non-electric surface friction and reducing mobility [21]. 
However, conflicting experimental results have been re¬ 
ported in which mobility is observed to increase with con¬ 
finement [20l[22l|23]. 

In order to resolve this discrepancy, we explicitly study 
the validity of the free-draining assumption for a freely- 
jointed charged chain that is squeezed by a radial po¬ 
tential that acts transverse to an electric field E = Ez 
(Figurej^. Our study does away with impermeable walls 
that add complications by introducing shear stresses, 
modifying friction coefficients, generating electro-osmotic 
flows and screening hydrodynamic interactions [24H26] . 
Only the effects of confinement on the polymer con¬ 
formations are preserved. Unless otherwise stated, we 
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FIG. 2. Electrophoretic mobility of confined polyelectrolytes 
normalized by the free-solution value /xo as a function of effec¬ 
tive confinement radius i^eff = y/k^T/k for various degrees of 
polymerization N. Inset (a): Asymmetry ratio (p = Rgz/Rgr 
as a function of i^eff normalized by contour length L. Open 
symbols denote neutral chains. Inset (b): Mobility as a func¬ 
tion of the asymmetry ratio p. 


use a harmonic potential 1/2 = kr‘^/2 , where k is the 
confinement strength and r is the distance from the z- 
axis. The ra dial pot ential forms a tube of effective ra¬ 
dius i^eff = y^kBT/k that imposes lateral constraints on 
the polymer conformation but does not act on the fluid. 
Although the potential confines the chain, no net axial 
force acts on the squeezed polyelectrolyte. Hence, the 
electro-hydrodynamic equivalence principle m Ell [28] 
might lead us to expect the drift velocity to match the 
free-draining value jHoE. 

Our simulations use the numerical methods reported 
by Hickey et al m- N purely repulsive Lennard-Jones 
charged beads (diameter a) are linked with finitely ex¬ 
tensible non-linear elastic bonds into a polymer of con¬ 
tour length L = (0.965cr) N. Monomers interact via a 
Debye-Hiickel potential and are embedded in a multi¬ 
particle collision dynamics (MPCD) fluid [30l |3T]. Un¬ 
less otherwise stated, the polymer chain is coupled to 
the MPCD fluid by including the monomers in each 
Andersen-MPCD collision event [32] [33]. The MPCD 
collision cells are cubic and of size a, which defines the 
unit of length. We employ periodic boundary conditions 
on a rectangular control volume of size 35 x 35 x 150 a^. 
MD monomers have a mass 3 times greater than the 
MPCD mass scale m and have one unit charge — le. 
The density of the fluid is 5/a^, while the monomer 
size is a = a/2, the Debye length is Ad = Icr and 
the Bjerum length is Ap = 1.5cr. The electric field is 
set to = 1 in simulation units [29], unless otherwise 
stated. Electro-hydrodynamic effects are implicitly sim¬ 
ulated by assigning MPCD particles a charge based on 
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EIG. 3. Scaled mobility p/po as a function of normalized 
effective conhnement radius Ref^/Rgo for various degrees of 
polymerization N and confining potentials Ui. Dotted lines 
represent the oriented-rod model and dashed lines the Long- 
Ajdari screened electro-hydrodynamics model {N =120 scar¬ 
let; A" = 80 amber; A = 40 beryl) The solid line shows Eq. [H 
Inset (a): Mobility dependence on {Res/L)‘^. Inset /Sj- 
Mobility dependence on (Ref^/Rgo)^, consistent with the form 
of Eq.[^ The fit produces an intercept of 3.0 ± 0.1 and slope 
of 15.3 ±0.3. 


the Debye-Hiickel approximation (Figure such that 
the charge density a distance r away from a monomer 
is Pe(^) cx: exp {—r/X]j)/r. If the resulting charge on any 
given fluid particle exceeds a threshold value, it is re¬ 
duced for both the MPCD and the associated monomers 
to account for charge condensation. This mean-field 
MPCD-MD Debye-Hiickel algorithm reproduces the non¬ 
monotonic increase in mobility with respect to A of 
charged oligomers, as well as the nonzero mobility of cer¬ 
tain net-neutral block polyampholytes [29] . 

The conformational changes due to confinement can 
be characterized by the asymmetry ratio p^ which is the 
ratio of the axial radius of gyration Rg^ to the radial 
component Rgr. When the harmonic potential is weak 
(Aeff larger than the unconfined radius of gyration Rgo)^ 
the chain is unperturbed and p = 1 (Figure [^). When 
the effective tube size is small, the polymer deforms into 
a string of de Gennes blobs [35]. Figurej^ shows that the 
asymmetry ratio of a polyelectrolyte radially confined by 
a harmonic potential scales as ^ L/Aeff, as theoreti¬ 
cally predicted for neutral chains [36] . 

For weak confinements, the mean-field MPCD-MD 
Debye-Hiickel simulations verify that the mobility is in¬ 
dependent of effective tube diameter and chain length 
remaining near the free-solution value /xq, as expected 
(Figure]^. However, in stronger confinements the elec¬ 
trophoretic mobility rises sharply above /xq. The increase 
in mobility above the free solution value appears roughly 
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FIG. 4. Velocity of the MPCD fluid about the polyelectrolyte’s centre of mass normalized by the electrophoretic velocity of 
the free chain. The degree of polymerization is iV = 200, the electric field is E' = 2 and the MPCD-in-MD coupling scheme 
is used [34]. Inset (a): Fluid velocity field for an unconfined polyelectrolyte. Inset (b): Velocity field for confinement by an 
axially symmetric harmonic potential with /c = 0.3 and an effective radius of Eeff = 1.8 (in units of MPCD collision cell size 
a). Inset (e): The far-held huid velocity held of Figure |^. 


inversely proportional to Eeff in Figure]^ The inset (Fig¬ 
ure If)’ however, reveals that this linearity is only ap¬ 
parently true for small asymmetries. Saturation to an 
asymptotic value of ~ 4/3 occurs at large (p. 

A semi-logarithmic representation of the elec¬ 
trophoretic mobility as a function of effective conhnement 
radius is presented in Figure Also included in Figure 
are simulations that replace U 2 with the steeper conhning 
potential Uq = kr^ /6 that forms a tube of effective size 
Eeff = ^k-sT/k and demonstrate that the mobility rise 
does not depend on the harmonic nature of the conhning 
potential. While Figuresuggests mobility dependence 
on asymmetry ratio, it obscures a length dependence for 
a given asymmetry ratio. Figure It demonstrates that 
moderate axial conhnement may lead to successful sep¬ 
aration of polylectrolytes for a given asymmetry ratio. 
However, actual experiments are done at a hxed value 
of Eeff, not at a hxed value of p. Normalization by the 
unconhned radius of gyration Rgo collapses the curves 
(Figure]^) and shows that mobility can be ht using an 
empirical function of the form 


Mo 
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We note that Eq.[^and Figure predict a plateau value 
li/po = 4/3 for strong conhnement. The dependence on 


the ratio {Ref^/Rgo)‘^ is not a result of the harmonic na¬ 
ture of U 2 since the proposed function (Eq. [B remains 
predictive for the Uq case as well (Eigure [^y Clearly, 
our data and Eq. support the idea that longer polylec- 
trolytes should migrate faster under conhnement; how¬ 
ever, this does not explain the physical mechanisms that 
leads to this potentially useful Wdependent rise in mo¬ 
bility. 

Length dependent electrophoretic mobilities often sug¬ 
gest that a chain is no longer free-draining [12]. How¬ 
ever, we hnd that this is not the case — the polyelec¬ 
trolyte remains free-draining. If it did not, long-range 
electro-hydrodynamic interactions would couple distant 
segments and the polyelectrolyte would electrophorese as 
a Zimmian chain, entraining distant huid along with it. 
We explicitly consider the entrainment of huid in Eig¬ 
ure In the absence of conhnement (Eigure [^), a small 
electro-osmotic how (magnitude < 5% of poE) occurs 
within the polymer coil. This is consistent with the ideal 
picture of a long free-draining polymer for which hydro- 
dynamic interactions are negligible and there are no far- 
held perturbations to the surrounding huid. 

The how-held predicted by our simulations of an elec- 
trophoresing chain conhned to an ehective tube of Eeff = 
1.8 is quite similar (Eigure [^). If free-draining were 
thwarted by conhnement then the entrained huid would 
be expected to move with a velocity comparable to the 
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Effective Confinement Radius, R^^/Rgo 


FIG. 5. The effective charge of confined polyelectrolytes 
(normalized by the free-solution value Qeff o) a function of 
normalized effective conhnement radius Ref^/Rgo for various 
degrees of polymerization N. Inset (a): Radius of gyration 
Rg of conhned chains. Normalization of Reff with RgQ ~ 
causes the confinement compactification minima to concur. 
Inset (h): The effective friction coefficient of confined chains 
normalized by the free-solution value ^eff o- Simulations with¬ 
out counterion condensation (hydrodynamic interactions) are 
labelled “no CIC” (“no HI”). Inset (c): The effective friction 
coefficient ^eff normalized by the rod-model ^rod, which shows 
the correction factor due to charge condensation. 


electrophoretic velocity of the chain and long-range per¬ 
turbations akin to those produced by a no-slip, imper¬ 
meable body would be expected. However, this is not 
observed. Rather, the entrained fluid velocity remains 
significantly less than that of the translating chain (Fig¬ 
ure IT)- Since neither the unconfined nor the confined 
coils produce significant far-field flows, they can both be 
described as “free-draining”. This is emphasized in Fig¬ 
ure Ef, which explicitly shows that far from the confined 
polyelectrolyte the velocity field is indeed zero. 

Since the electrophoresing chain remains free-draining, 
the mobility must be a function of effective local terms. 
Because electrostatic repulsion is not fully screened 
within the Debye layer, local segments are locally stiff 
and form roughly rigid rods of length ^ Ad- In free- 
solution, these local segments randomly sample all ori¬ 
entations but when the chain is strongly confined the 
segments orient and the effective friction coefficient of 
the segments decreases to approximately ^eff ~ ^rod ~ 
^11 (5 — 2S) /3, where S = (3cos^ 0 — l) {2 is the ori¬ 
entation order parameter and is the friction coeffi¬ 
cient of a slender rod oriented parallel to z. This sim¬ 
ple, local rod-model predicts that the electrophoretic 
mobility should increase as ^ which agrees qual¬ 

itatively with MPCD-MD simulations although it does 
not predict any A/'-dependence exhibited in simulations 


(Figure dotted lines). Likewise, more meticulously 
modelling of monomers as screened point-charges by the 
Long-Ajdari electro-hydrodynamic interaction tensor [13] 
produces more accurate predictions but the mobility re¬ 
mains A^-independent (Figure]^ dashed lines). 

This leads to the conclusion that the A^-dependence 
must be a result of variation of the effective charge Qeff, 
which occurs as a higher order effect to the primary 
mechanism of A/"-independent friction reduction due to 
orientation along the tube. To test this explicitly, the 
counterion condensation is “switched off” [29| such that 
the mobility is determined solely by the effective friction 
of a segment. Measuring the effective friction directly 
from the mobility, we observe that ^eff is well approxi¬ 
mated by ^rod (Figure [^; no CIC). Thus, local friction 
coefficient reduction due to orientation within the Debye 
layer is indeed the main mechanism for the mobility in¬ 
crease; however, it is not sufficient to account for the N- 
dependence. If, instead, hydrodynamic interactions are 
“switched off” |37l |38] then the mobility still varies from 
its free-draining value. Since the friction coefficient re¬ 
mains constant by construction, the effective charge must 
vary. Figure]^ (no HI) shows the explicitly measured Qeff 
as a function of confinement radius. An A/'-dependent 
minimum of Qeff exists when hydrodynamic interactions 
are included (Figure [^. It is this A^-dependence in Qeff 
that leads to the A/'-dependent mobility. 

The effective charge varies due to confinement-induced 
variation in monomer density and, in particular, to the 
A/'-dependent phenomenon of confinement compactifica- 
tion [39]. At moderate confinements, the chain first ori¬ 
ents along ^ causing Rg to decrease (Figure]^). Since the 
likelihood of finding distant segments near one another 
increases, diffuse layers overlap (Figure Q and counte¬ 
rion condensation increases, which lowers the effective 
charge relative to the unconfined coil (Figure]^. More 
strongly confined chains are stretched into a string of 
blobs along the axis of the tube causing Rg to increase 
(Figure [^). Blob theory predicts that the density of 
monomers continues to increase with confinement. How¬ 
ever, in the strong confinement limit, the chain is highly 
stretched and no longer in the blob regime, such that dis¬ 
tant segments are then less likely to be near one another. 
In turn, the fraction of condensed counterions decreases 
and Qeff increases toward the charged-rod limit m- 
The position of the Rg minimum universally depends on 
R^e/RgO (Figure!^) [H] and, therefore, so 

too does Qeff (^eff/^^o)- II is precisely this secondary 
effect that causes electrophoretic mobility /i to be a uni¬ 
versal function of the effective tube radius when normal¬ 
ized by unconfined radius of gyration {Ref^/Rgo) as seen 
in Figure 

Finally, counterion condensation has a tertiary effect 
on /i. When more counterions condense increased cou¬ 
pling occurs and the effective friction coefficient decreases 
since the effective charge sets the number of counterions 
remaining in the diffuse layer screening hydrodynamic 
coupling. Likewise, when the effective charge goes up, 
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hydrodynamic screening is enhanced and ^eff rises (Fig¬ 
ure [^). This effect ensures that the mobility is not quite 
Qeff/Crod since a correction factor c((5eff) = Ceff/^rod 
must be accounted for in the effective drag (Figure]^). 
Thus, we conclude that the mobility increases from its 
free-draining value primarily because orientation reduces 
the local friction coefficient, while length dependence 
arises from the secondary effect of charge condensation 
due to confinement compactification, which in turn cause 
the tertiary effect of reduced hydrodynamic screening, 
further decreasing the effective friction coefficient. 

By considering a simplified system that confines freely- 
jointed polyelectrolytes via radial potentials rather than 
impermeable walls, this study explicitly demonstrates 
that the electrophoretic mobility depends on conforma¬ 
tion in a length-dependent manner, though the chain, 
as a whole, remains free-draining. The primary effect is 
local alignment of segments within Debye layers, which 
would not be accounted for by blob theories of confined 
polymers that assume uniform monomer density and ori¬ 


entation distributions within blobs [42]. Simulations of 
this effect require computational techniques that account 
for finite Debye layers, such as the mean-field MPCD-MD 
Debye-Hiickel algorithm used here [29] . 

Our results demonstrate that in microfiuidic devices, 
frictional drag with walls is not entirely responsible 
for changing mobility — wall-chain interactions increase 
drag competing with alignment and counterion conden¬ 
sation effects. We expect that our fundamental findings 
will aid the design of new electrophoretic methods for 
highly charged biomolecules. While the polyelectrolytes 
in this work are freely-jointed with a finite Debye length, 
DNA has a large persistence length. Radial confinement 
of DNA would still orient local segments; however, N- 
dependent charge condensation due to overlapping Debye 
layers is predicted to be less significant. 
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